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Abstract 

The perturbative construction of the 5-matrix in the causal spacetime 
approach of Epstein and Glaser may be interpreted as a method of regu- 
larization for divergent Feynman diagrams. The results of any method of 
regularization must be equivalent to those obtained from the Epstein- Glaser 
(EG) construction, within the freedom left by the latter. In particular, the 
conceptually well-defined approach of Bogoliubov, Parasiuk, Hepp, and Zim- 
mermann (BPHZ), though conceptually different from EG, meets this re- 
quirement. Based on this equivalence we propose a modified BPHZ procedure 
which provides a significant simplification of the techniques of perturbation 
theory, and which applies equally well to standard quantum field theory and 
to chiral theories. We illustrate the proposed method by a number of ex- 
amples of various orders in perturbation theory. At the level of multi-loop 
diagrams we confirm that subdiagrams as classified by Zimmermann's forest 
formula in BPHZ can be restricted to subdiagrams in the sense of Epstein- 
Glaser, thus entailing an important reduction of actual computations. The 
relationship of our approach to the method of dimensional regularization 
(and renormalization) is particularly transparent, without having to invoke 
analytic continuation to unphysical spacetime dimension. It sheds new light 
on the role of some parameters that appear within dimensional regulariza- 
tion, and thus establishes a direct link of this traditional method to the 
BPHZ scheme. 
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1 Introduction 



As is well-known, formal perturbation theory applied to relativistic quantum field 
theory, in general, leads to ill-defined expressions for elements of the ^-Matrix. 
Integrals that are expected to describe probability amplitudes for certain scat- 
tering processes are found to be divergent. Generally speaking, these ultra-violet 
(UV) divergencies can be traced back to the naive application of time-ordering 
in the description of propagation of particles. There are two conceptually rather 
different lines of attack to deal with this problem: The first of these consists in a 
set of regularization procedures all of which are designed such as to replace diver- 
gent integrals in Feynman diagrams by convergent ones in a consistent manner. 
These empirical regularization schemes are justified by their usefulness in practical 
applications of quantum field theory to physical processes. In order to be consis- 
tent, they must fulfill all physical (normalization) conditions, order by order, or, 
at least, must contain enough freedom to meet these conditions. This is the es- 
sential prerequisite for the procedure of renormalization. In other terms, not every 
scheme of regularization of divergent integrals of a given theory meets the stronger 
requirement of renormalizability of that theory. 

The second line follows the approach developed by H. Epstein and V. Glaser, [1], 
which is based on causality and locality in coordinate space. This procedure makes 
use of a well-defined rule for time ordering and thereby allows to construct an 
entirely divergence-free (S-matrix from basic and general principles. The Epstein- 
Glaser (EG) approach is mathematically rigorous, within perturbation theory, but, 
when applied without modifications, is not very useful in practice. By its very 
construction, due to the process of distribution splitting, it contains a certain 
freedom which, subsequently, is fixed through its interpretation in terms of physics 
in the process of renormalization. In fact, as was first proposed in [2, 3], the EG 
method can be interpreted itself as a regularization scheme. Thus, this approach is 
particularly useful as a reference framework for testing whether a given empirical 
method of regularization is physically admissible (in the sense of renormalizability) , 
or not. 

Among the regularization procedures of the first group the classical method of 
Bogoliubov, Parasiuk, Hepp and Zimmermann (BPHZ), by its conceptual clarity, 
is the most rigorous one [4, 5]. The rules of regularization that it contains at the 
level of one-loop diagrams are equivalent to those of EG, with the latter suitably 
translated to momentum space. In addition BPHZ contains a general prescription, 
in the form of the forest formula, for regularizing higher-loop diagrams. However, 
although its logical structure is transparent, the BPHZ procedure leads to rather 
involved integrals in explicit calculations which make it less suitable for practical 
computations of Feynman amplitudes as compared to more empirical methods 
such as dimensional regularization or the like. (In using the term "dimensional 
regularization" we follow common conventions. In fact, this nomenclature means 
dimensional renormalization with, say, minimal subtraction.) 
With d{x) a scalar distribution of singular order u the EG method defines ad- 
vanced and retarded distributions through splitting of its support by a space-like 



2 



hypersurface, say v ■ x — with v a timelike vector. As is well-known, a construc- 
tion valid for all physically relevant values cu and for all test functions g e ©'(M'^), 
is then 

J d^x d^et,reg{x)g{x) = J d'^x d{x)Q{v ■ x) {Wg) (x) , (la) 

J d^X dadv,reg(2:)^(x) = " /d^^ d(x) [1 - e{v ■ x)] (W g) (x) , (lb) 

where the operator W is defined through its action on test functions g{x) 

^ a 

{Wg) {x) = g{x) - w{x) ^ (^"^) (2a) 

Qj. 

\a\=0 



with the customary short-hand for partial derivatives 

Qai+...+ak 

dxl' ■ ■ ■ dxl" 



D'-TT^, \a\^ai + ... + ak, 



with the standard multicomponent notation for the coordinates, and with w{x) 
a function satisfying the conditions 

w(0) = 1 , D'^wl^^Q = for all 1 < |a| < . (2b) 

A given choice of the function w represents a specific regularization. This is the 
perspective adopted in [2, 3]. Any two different regularizations differ by a S- 
distribution and derivatives thereof up to the order lv in the integrands, [6], viz. 

d'^x {d^eg,wi{x) - d^eg,w2{x)) 9ix) = (i^x ^ CaD^dix) g{x) with 

\|a|=0 / 



I X 



Ca = jd^'x d{x){-lt\— [w2{x) - W^{X)] . 

This freedom of choice is essential for the subsequent renormalization process which 
relates the free parameters to values of physical parameters. As also shown in 
that work, a modified subtraction operator such as the one proposed in [7] will 
yield a valid regularization but may turn out to be too restrictive for successful 
renormalization. 

In this paper we study a new method that we propose to call modified BPHZ proce- 
dure. This method combines the practical usefulness of dimensional regularization 
with the structural simplicity of the classical BPHZ renormalization in the light 
of its equivalence to the Epstein-Glaser method. Like for any other empirical reg- 
ularization method, the rigorous Epstein-Glaser framework is the landmark with 
respect to which the correctness and use of our modified procedure must be judged. 
The paper is organized as follows. In Sect. 2 we discuss the equivalence between 
the BPHZ and EG frameworks. In Sect. 3 we describe the idea of the modified 
BPHZ method and its justification, by means of its relationship to Epstein-Glaser. 
In Sect. 4 we give some instructive examples and work out the relationship to 
dimensional regularization. The final Sect. 5 gives a summary and outlook. 
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2 Equivalence of BPHZ and of EG frameworks 



The BPHZ scheme is based on Feynman rules in momentum space. Schematically, 
and at this point still somewhat formally, a given diagram 7 with internal momenta 
k is translated to an integrand of the form 



(3) 



lec 



vev 



where p denotes the set of external momenta, while k stands for the internal 
momenta to be integrated over. The factors Ac are proportional to Feynman pro- 
pagators Ai;' in momentum space, and correspond to the internal lines / of a given 
set jC. The momentum flow is defined by the conventions chosen in the forest 
formula. A given vertex V in the set V of vertices contributes the factor Py. 
Consider an arbitrary irreducible one-loop diagram whose degree of divergence is 
(i(7). The BPHZ approach replaces the integrand by the modified expression 



R,ip,k)^{l-4^^^)l,ip,k) , 
where S''^ = V - f 



d 



|n|=0 



dp" 



(4a) 
(4b) 



The Taylor operator stands symbolically for the expansion in terms of the set 
of independent external momenta p. In view of subsequent renormalization, the 
general result of regularization has the form 



d^kl^{p,k) 



BPHZ,reg 



jd^kR^{p,k) + P^'^^''^\p) , 



(5) 



with p(°'(T))(p) a polynomial of degree ci(7) representing the remaining freedom. 
For the sake of illustration we will refer repeatedly to 0^ theory in which case Py 
yields a power of the coupling constant g. As an example, consider the four-point 
function of 0^ theory at one loop i.e. the diagram shown in Figure 1, with external 
momentum p and containing two internal lines. In this case (i(7) — and BPHZ 
regularization yields 




Figure 1: Four-point function in the 0^ model, at one loop, with external momen- 
tum p = Pl + P2 



A|(P) 



BPHZ.reg 



(2vr) 



d^A; 



(6) 
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In order to compare with the corresponding result of EG regularization the equa- 
tions (la) and (2a), as well as all test functions are transformed to momentum 
space, by Fourier transform, so as to obtain 

dregig) = J d^a; dreg{x)g{x) = J d'^k dreg{k)g{k) , (7) 



as well as the analogue of (6), 

reg (27r)6 J k^ ~ {{p-ky- 



' W ... T } ■ (8) 



In the present example one may choose the function w tohew{x) = 1, as a limiting 
case. Obviously, it satisfies the conditions (2b). Furthermore, its Fourier transform 
being w{p) — l/{2TryS{p), it is seen to yield the BPHZ expression (6). 
Of course, proving the equivalence of the BPHZ and EG schemes, beyond the 
one-loop level and for other theories, becomes technically more complicated. As is 
well-known, this is due to the fact that EG is an expansion in terms of the number 
n of vertices, i.e. in terms of powers of the coupling constant, while BPHZ is an 
expansion in terms of loops, i.e. a formal expansion in terms of Planck's constant. 



3 A modified BPHZ procedure 

In the classical BPHZ method divergent momentum integrals are regularized by 
means of appropriate Taylor subtractions of the integrand. Although, on the basis 
of Zimmermann's forest formula, this approach is transparent and well-defined in 
principle, its practical implementation in higher orders of perturbation theory is 
cumbersome. From a practical point of view, other methods of regularization such 
as analytic continuation in the dimension of space-time are better tools in actual 
calculations. 

The alternative procedure that we propose aims at modifying the well-defined 
framework of BPHZ in such a way that it becomes as practicable as customary di- 
mensional regularization. In essence, the idea is to introduce Feynman parameters 
at the level of the unsubtracted integrand and to apply Taylor subtraction to the 
modified integrand only. As we will show this preserves the mathematical rigour of 
the BPHZ scheme but simplifies enormously subsequent integrations over internal 
momenta. 

To start with and in order to explain the essence of the modified method we give a 
very simple example from the (f)^ model. Within the BPHZ scheme and at second 
order in the coupling constant g, the contribution of the one-loop diagram to 
the four-point function is logarithmically divergent. BPHZ regularize it by Taylor 
subtraction of the integrand to order zero, viz. 
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and thus obtain a well-defined expression. In a first step, and in analogy to di- 
mensional regularization, we parametrize the unmodified integrand by means of 
a Feynman parameter z. such that the integrations over z and over the internal 
momentum A; may be interchanged. In a second step the integration variable is 
subject to a translation by the vector [z — k'r^q = k\{z — l)p, so that 
mixed terms containing external and internal momenta no longer appear. Finally, 
the Taylor subtraction is applied to the modified integrand. The three steps are 
given by, respectively, 

m = T^-n [d'k (1 - ti) [dz - — — — ^ — -2 

^ 2(lo^ L^'j"^'^ " [k^-2pk{l-z)+p\l-z)- 



Note that the scaling behaviour of the integrand for large values of k remains 
unchanged by the introduction of a Feynman parameter. Therefore, the coefficients 
of the Taylor expansion of order higher than the degree of divergency lead to 
convergent integrals. The calculation of the integral over q is standard. Making 
use of a Wick rotation, one obtains 

Note that, unlike in dimensional regularization, the result (11a) is exclusively ob- 
tained in dimension four of physical space-time. Furthermore, the remaining free- 
dom in the approach discussed here which allows for a constant additive term (with 
respect to p), may be made explicit by replacing m ji, with jj, an arbitrary mass, 
in the numerator of the logarithm in (11a), i.e. A(p) may be replaced by 

Indeed, the expressions (11a) and (lib) differ by a constant only. For instance, 
the specific choice //^ = ^^'dira.xeg^~'^ ^ with 7 Euler's constant, reproduces the 
well-known result of dimensional regularization. 

Somewhat more generally, a convergent one-loop integral whose integrand was 
Taylor subtracted to the appropriate order u) 

J,{p)^ Jd'k{l-t^)l^{k,p), (12a) 
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is transformed by a translation of the argument, k t-^ q — k + Xp, 

J,{p)^ Jd'k {{l-t^)l^{k + Xp,p)-A{q,p)} , (12b) 

where the function A{q,p) is defined by this equation and, hence, is given by 

A{q = k + Xp,p) = {l-t;:)l,{k + Xp,p)- . (13a) 

Note that in the first term the Taylor operator applies to both arguments of the 
function while in the second term it applies to the second argument only. 
Alternatively, the function A(g,p), with q = k + Xp can also be written as follows: 

A{q = k + Xp,p) = {l-t-) [(t-A(g,p))|^^,^J . (13b) 

The result (13b) follows from an identity for the Taylor operator tp in the variable 
p about the point p — 0, applied to a diff'erentiable function F of two variables, 

t-pF{k + Xp,p) = (t,-F(g,p)|,.,^,,) . (14) 

Indeed, denoting by di and 82 the derivatives with respect to the first and second 
argument of F, respectively, the left-hand side is 

t^Fix + Xy,y)^J2{N) {9idt'F{x, 0)) . 

i=0 ^ ^ 

The right-hand side, in turn, is computed to be 

f ^ 1 



, fc=0 



N / . X N . 

\ ^ / t \ \ 1 



With {d'^~^y^)\ = k\5N-ik this is seen to be the same expression as above. 
As a result, the integral J-y, after translation of the internal momentum, takes the 
form 

J,{p) = Jd'k (1 - t^) I,{k + Xp,p) + jW(p) , (15) 

where J^^^ is the integral over A. Closer examination of (13a) shows that this 
integrand can be written as a sum of derivatives with respect to k of order one and 
higher, and, hence, gives rise to surface terms which vanish at infinity. Thus, J^^^ 
vanishes. This calculation demonstrates that translation of the integration variable 
is an admissible operation. 

As will be clear from the examples worked out below, the parameter A, in gen- 
eral, is a function of the Feynman parameter(s) z. Like in the example above the 
translation is chosen such that mixed terms in k and p disappear. The integrand is 
then a function of k'^ only so that the integration can be done in Euclidean polar 
coordinates, via Wick rotation. 

These examples motivate the following modified BPHZ procedure: 
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1. In a given integral Ij{k,p) with external and internal momenta p and k, 
respectively, introduce integral representations by means of a set of Feynman 
parameters z, interchange the z integrations with the operator (1 — tp), and 
with the integration over the internal momenta k. 

2. Perform a translation of the A;- variables such that internal and external mo- 
menta are decoupled. 

3. To the integrand apply Taylor subtraction up to singular order cu. 

4. Do the fc-integrals by means of Wick rotation and using Euclidean polar 
coordinates. 

5. In order to make contact with dimensional regularization replace the mass 
parameter (s) by general constant (s) n so that only additive terms appear 
which form a polynomial in p up to and including the singular order u. In 
some cases this does not exhaust the freedom necessary for renormalization 
because, obviously, the modified BPHZ method has the same number of 
parameters as the original one. This is essential for identifying the physical 
parameters of the theory (masses, charges etc) in each order of perturbation 
theory. 

This is a well-defined algorithm whose advantages are evident. The general analysis 
given in eqs. (12a)-(15) as well as the examples at second and higher orders, 
lend strong support to the conjecture that it meets all requirements of physical 
renormalization. Its closeness to the original BPHZ regularization and, hence, to 
Epstein-Glaser regularization, guarantees that it is an admissible regularization 
scheme. 

By a suitable choice of the parameter(s) n one makes contact with well-known 
regularization methods such as dimensional regularization, without having to con- 
tinue to unphysical space-time dimensions. The method is mathematically rigorous 
but more practicable than the original BPHZ approach. 

Furthermore, turning to fermions, no continuation of the Clifford algebra of Dirac 
7-matrices is necessary given the fact that the modified BPHZ method works 
exclusively in dimension four. 

Our approach is rather close to the framework of Epstein and Glaser, but allows for 
direct comparison with unmodified BPHZ. Due to cancellations of a certain class of 
subdiagrams there are important simplifications in the calculation of higher-order 
processes. In order to explain this point we need some preparation and definitions. 
As we stated above, EG is an expansion in terms of the coupling constant g, while 
BPHZ is an expansion in powers of h, hence in terms of the number of loops. 
EG constructs a functional T„ describing a diagram with n vertices by recurrence 
from the tempered distribution Ti = i£int- The total diagram depends on fimc- 
tionals which were regularized previously at orders lower than n, say m < n. Thus, 
the corresponding subdiagrams contain irreducible parts with a number of vertices 
smaller than n. We shall call such subdiagrams Epstein- Glaser-suhdiagrams or, for 
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short, EG subdiagrams. The BPHZ framework, in turn, works by successive addi- 
tion of counter terms proportional to ascending powers of h and, as a consequence, 
requires a different classification of subdiagrams. Let us call BPHZ subdiagram any 
irreducible divergent part of the total diagram which contains a smaller number of 
loops than the main diagram. In particular, there will be subdiagrams which are 
lower in loop order but do not have a smaller number of vertices. We call these pure 
BPHZ subdiagrams. An example we shall study in more detail below is the "sun- 
rise" diagram in the model, cf. Figure 2. In the framework of BPHZ it contains 
three logarithmically divergent subdiagrams. In the perspective of EG, in contrast, 
it is a diagram with two vertices and, hence, contains no divergent subdiagram at 
all. In our terminology the three BPHZ subdiagrams are pure BPHZ subdiagrams. 
The sum of the counter terms generated by these subdiagrams does not contribute 
to the regularization of the sunrise diagram. This example as well as other exam- 
ples studied in [6] confirm this to be a general rule, and are in accordance with 
a theorem by Zimmermann [10]: Pure BPHZ subdiagrams do not yield counter 
terms, i.e. their sum vanishes, and, thus, they may be left out in the modified 
approach. 



p-k-q 




Figure 2: Sunrise diagram in the 0^ model 



We illustrate the method by a number of significant examples in second and higher 
orders. 



4 Examples 

We start with some classical examples from quantum electrodynamics and elec- 
troweak interactions, the self-energy of the electron, the vacuum polarization, and 
the vertex correction at one-loop order, then mention briefiy the case of the triangle 
anomaly. We finish with a typical second-order, two-loop, process and with some 
remarks about higher-order processes which illustrate the simplicity of our alter- 
native scheme. In all these examples the equivalence to EG regularization proves 
the correctness of the modified BPHZ approach. 

4.1 Quantum electrodynamics with electrons 

In the original BPHZ framework the self- energy of the electron reads 
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In the modified BPHZ approach we introduce a Feynman parameter interchange 
integrations, and substitute k ^ q — k — zp such as to decouple internal and 
external momenta, to obtain 

This is easily worked out to be 

The remaining freedom is made explicit by replacing by an arbitrary squared 
mass //^, 

One verifies that the choice 

/^^ = 47r/.^i..ege(^M (17) 

reproduces (the finite part of) the result known from dimensional regularization, 
see e.g. [11]. 

Lowest order vacuum polarization in original BPHZ is given by 

n..W = ji\ (1 - '3 (7>#3^T. (;_;,_^, j (18a) 

In the modified scheme, we introduce a Feynman parameter z, interchange the 
integration over z with the one over the internal momentum q, and perform a 
translation of the integration variable q^q — q — — z)p, to obtain 

n!.^W = -^/d. (s,./-p„p„).(l-.)ln(-,-j^!-^) . (18b) 

As before, in order to exhaust the remaining freedom, we have replaced the nu- 
merator in the logarithm by an arbitrary squared mass /i^. The (finite part of) 
the known result of dimensional regularization [11] is recovered by the choice 

/^^ = 47r/^dim.reg (19) 

The vertex correction, at the same order, finally, is found to be 

167r2 /^^ / ^ ''■'^''^^"^ ~ - 2;^ + m]7«[(l - x)^ - yf + m]7''} , 

(20) 
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where the denominator in the integrands stands for 



D — [x + y)m — x[l — x)p — y{l — y)p + 2pp xy . 

As before, we replaced the numerator in the logarithm by an arbitrary term 
/x^, to cope with the remaining freedom after regularization. The analogous re- 
sult in dimensional regularization is recovered by the same choice (17) as for the 
self-energy. This shows that the modified BPHZ regularization fulfills the Ward- 
Takahashi identity 

^E(p) = -A„(p,p). (21) 
Though not surprising, this is a consistency check. 



4.2 Chiral anomaly 



We also analyzed the well-known vector-vector-axial vector (VVA) chiral anomaly 
sketched in Figure 3, within the modified BPHZ procedure. 




Figure 3: Triangle graph contributing to the anomaly 



Denoting the amplitude by Tq^j, and choosing the internal loop momenta as shown 
in Figure 3, conservation of the vector current at the two lower vertices should yield 
the Ward identities 

p''T^^.u - , q'T^^, = , (22a) 

whereas the axial current vertex should produce an anomalous Ward identity which 
survives even in the limit of the fermion mass m going to zero, viz. 

(p + qTTa^^u = 2mT^, + ^e^.^r^V , (22b) 
the term T^,^ being given by 

^^'^^IWr lfc^^^ (fc-p)^-m^ ^ ^fc + g)^-m^ ^- j +(P - ^, /. - .) . 

11 



It is known that the anomaly can be shifted from the axial vector current to the 
vector current, or to a linear combination of these [8], [9]. Thus, by requiring that 
it be the vector current which is conserved, some of the freedom in the renormal- 
ization process is made use of. 

The diagram of Figure 3 is linearly divergent. If regularized by Taylor subtraction, 
in the spirit of BPHZ, it is given by 

IT ^ _ /"^ n _ ^1 )tr r 7 f-i^ + m ^ + i + m 

tr [{^ + m)7^(^: - + m)'-ial-,{f- + i + m)7^] 

[[k^ — m?){\ — X — y) + {{k — pY — m?-)x + {{k + qY — m?)yf' 

1 .1-. . ^4^ 



JQ 



-^Idxl dy/^(l-4) 



tr [{jt + m)-fu,{^ - ^ + m)jal5{jt + i + m)-fu] 



[{k + {qy — px)Y — {qy — pxY ~ m?{^ — x — y) -\- {q^ — m?)y + (p^ — m?)x\ 

A second term contributing to the chiral anomaly is obtained by interchanging 
(p g), (// i/). By the symmetry of the integrands this second term yields the 
same result as the first so that the factor 1/2 on the left-hand side can be dropped. 
Following the rules of the modified BPHZ scheme described in Sect. 3 one performs 
the substitution 

k — k — {qy — px) 

such as to separate internal and external momenta, and to allow for separation of 
terms even and odd in the new integration variable k. Indeed, only the even terms 
contribute to the integral. A straightforward calculation leads to the result 

the logarithmically divergent term and the finite term being given by, respectively, 
= ^ IdxJ^' \ye^,,, {(3a; - 1^ - {^V ' 1)?^ 



In 



2 

m 



m? + {qy — pxY — q^y — p^x 



+ 



+ 2 



27r2 7o ^0 

egf^va {[{y - l)q'' -{x- l)p''] [{qy - pxf - m^] - yq^p" + xp^q"} 
m? + {qy — pxY — q^y — p^x 

ySanarP^q^'Ky - 1)^1. - Xp^] + XeauaTq''p''[{x - l)p^ - 



+ {qy — pxY — Q^y — P^x 
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In this example the well-known result obtained from unmodified BPHZ, or from 
dimensional regularization [9], is obtained in a straightforward and technically sim- 
pler fashion. Replacing in the logarithmic integrand by an arbitrary parameter 
does not change the total expression. Indeed, this replacement produces an 
additive term proportional to 



{(3x-lK-(3y-l)?'^}ln(^^) , 



which yields zero after integrating over the Feynman parameter y from to (1 — a;), 
and over x from to 1. 

For the sake of completeness we verify the Ward identities (22a) and calculate the 
anomaly (22b). Straightforward calculation of the divergence p^Taij,v leads to the 
result 

P Tctai/ = — :j£afii'aP Q I dx I dy ^ ■ -. -^^ ^ 

ivr^ Jo Jo + \Q.y ~ V^') ~ hU ~ F ^ 

\{—yx + 2yx^ + — \x'^)p^ + {yx^ — y^x)qp + (— y"^ — 2y^x + ^y^ + yx)g^] . 

In the diagram of Figure 3 the vector bosons at the lower vertices are identical so 
that p^ = (and equal to zero in the case of external photons). With (f = p^ 
the integrand is antisymmetric under exchange of x and y while the domain of 
integration is symmetric. Therefore, the integral vanishes and the first of the Ward 
identities (22a) holds true. The second Ward identity follows from the first by the 
symmetry {p ^ q), {fi ^ v). 

Regarding the divergence (22b) which contains the anomaly, we find for the first 
term on the right-hand side 

1 /"^"^ 

2mT^^{m) = i^e fj^varq" p'' dx dy — . 

TT^ Jo Jo + [Qy - PX) - TV - P^ 

The left-hand side of (22b) is calculated along the lines of the procedure described 
above. We find the result 

"1 /"l-a; , ^2 



and, upon comparison with the previous formula, 

(p + q^Ta^u = 2mT^^{m) + -^e f^y^rq" p'' , 
which is, indeed, the anomalous identity (22b). 

This example illustrates well the advantage of the modified BPHZ procedure, as 
compared to original BPHZ renormalization or to dimensional renormalization, 
by its simplification of the momentum integral. Furthermore, the equivalence to 
Epstein-Glaser regularization puts the modified procedure on solid ground. As 
compared to dimensional regularization, in particular, there is no need to introduce 
an analytic continuation of 75 to any other space-time dimension than four. 
We note in passing that the chiral limit m — > 0, like in the usual BPHZ framework, 
requires a separate discussion. We do not treat this case in the present work. 
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4.3 Higher loop diagrams 

The "sunrise" diagram in the 0^ model, cf. Fig. 2, provides an instructive example 
for the comparison of BPHZ and EG regularizations. Being a diagram with two 
vertices it contains no EG subdiagrams at all. In the perspective of BPHZ, however, 
it contains three logarithmically divergent pure BPHZ subdiagrams. Thus, in the 
former case it is regularized in a single step by Taylor subtraction with respect 
to the external momentum, while in the latter, one would have to invoke the 
forest formula for identifying the counter terms stemming from the three divergent 
subdiagrams. This is to say that the modified approach which, in essence, is a 
practicable version of EG, is technically simpler, and, furthermore it uses the fact 
that, in accord with Zimmermann's theorem [10], the contributions from all pure 
BPHZ subprocesses cancel. 

Regularizing the quadratically divergent diagram of Fig. 2 by Taylor subtraction 
of the integrand one has 

^^^^ " 6(lo« j'^'^^j^'^ - ^9 [(p_fc-g)2-^2] -grz—i -^fzr^2 ■ (23a) 

One succesively introduces Feynman parameters for the internal momenta /c, g, 
and p — k — q, and applies the necessary translations which decouple internal and 
external momenta. Details of this calculation are given in the appendix. The result 
is 

9' [' (l-2.)(l-2x)p^ 



{z-l){l-z + z- 



2) 



[zjl - z)il - x) + x]m^ \ ^ 

\ —xz{l — ^) (1 — x)p'^ + z{l — z){l — x)m? + xw? J 

As in the previous examples one replaces the parameter by an arbitrary pa- 
rameter fj? but verifies that the result (23b) remains unchanged, 

E('^)(p)=E(p). 

It is instructive to compare the result (23b) to a calculation of the sunrise diagram 
using dimensional regularization [12]. The result is 

= ^ j^iz ^d. { (-^3m^ + (1 - .)p^) 

[z{l — z){l — x) + x]m'^ \ 1 2 



\— a;z(l — z){l — x)p'^ + z{l — z){l — x)'w? + xm?) ~'~ 2^ J ^ ^ 

The expressions (23b) and (23c) both are regularizations of the same scalar dis- 
tribution. Furthermore, their Taylor expansion vanishes up to the order p^, in the 
first case by construction, in the second case due to the additional term p^/2. As 
this exhausts the remaining freedom in regularizing, one concludes that the two 
results are identical. 
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In order to make contact with the unmodified BPHZ procedure we have checked 
by explicit calculation using the forest formula that, indeed, the three subdiagrams 
which are not EG subdiagrams cancel in the Taylor expansion and, hence, do not 
contribute to the regularization of the sunrise diagram. Thus, the modified BPHZ 
procedure is very close to pure Epstein-Glaser regularization and avoids from the 
start irrelevant contributions from pure BPHZ subdiagrams to the regularized 
amplitudes, in agreement with the proof by Zimmermann [10]. 
These results are corroborated by case studies of EG regularization in higher orders. 
Among others we studied the four-point fimction of the 0^ theory in dimension 4, 
at the level of two loops. The same model in dimension 6 provides an example 
which besides yielding divergent EG subdiagrams, also exhibits pure BPHZ subdi- 
agrams. The contribution of the pure BPHZ subdiagrams, i.e. those which have no 
counterpart in EG, upon Taylor subtraction, are found to vanish, as expected. In 
the latter example we also studied three-loop contributions to the two-point fimc- 
tion [6]. In all cases EG regularization on one hand, and calculation following the 
forest formula restricted to EG subdiagrams, on the other, yield identical results. 

5 Conclusions and outlook 

The modified BPHZ procedure that we advocate in this paper combines the trans- 
parent concept of BPHZ regularization with the practical usefulness of dimensional 
regularization. In particular, the relevant integrals are easier to calculate than the 
corresponding ones within the original BPHZ method. Furthermore, the momen- 
tum dependent logarithms always contain a reference mass which is identical with 
the typical mass parameter of the theory (the electron mass in the case of the 
examples from QED, the scalar mass in the 0^ model). We showed, however, that 
rescaling is possible within the freedom allowed by the regularization process. In 
the examples with one loop, for instance, this allows to introduce a new mass pa- 
rameter which may be identified with the parameter of dimensional regularization. 
However, there is an essential difference here: While dimensional regularization 
requires the introduction of this parameter for (spacetime-)dimensional reasons, 
in our approach it is a manisfestation of the general freedom within the process 
of regularization. This remark, in turn, justifies its appearence in the results of 
dimensional regularization. 

The comparison of BPHZ regularization along the forest formula with the Epstein- 
Glaser construction confirms the expected significant simplification of explicit cal- 
culations in higher orders. In the light of different classifications of subdiagrams 
in the framework of BPHZ on one side, and in the Epstein-Glaser construction on 
the other, the summation over the subdiagrams contained in the forest formula is 
restricted to subdiagrams in the sense of Epstein-Glaser.^ The modified procedure 
implies, in particular, that the combinatorics of higher-order diagrams is described 
by the restricted forest formula which takes account exclusively of the class of 

^Of courso, a cortain choico of the standard momentum flow in the forest formula had to be 
made but the conclusion should be independent of that choice. 
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EG subdiagrams. Thus, this method is as straightforward as, e.g. dimensional 
ularization, and has the virtue to rest on solid mathematical ground. 
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A Derivation of Eq. (23b) 

The strategy for deriving (23b) goes as follows. A first Feynman parameter denoted 
by z is introduced for the /c-integration. A subsequent translation by k \^ k = 
k — {1 — z){p — q) then frees this inner momentum from mixed terms. Furthermore, 
we introduce the redundant operation (1 — t^), viz. 

m ■■= ^s(p) = jdu (1 - 4) j\\z Jd'k (1 - ti) 

1 1 



{(1 — z)[{p — k — qY — w?] + z{k'^ — m^)}^ q^ — m? 

1 



d^g (l-^p)/d. jd'~k {l-tD 



/pi POC 
(1 - ti) / dW dp (1 - ti) 
Jo Jo 



+ z{l -z){p- qf - m2}2 q^ - m? 
1 roo 3 1 



{p^ — z{l — z){p — qY + w?}"^ q^ — rri 



2 



Doing a partial integration in the integral over the parameter z, and introducing 
the abbreviation :— m^/{z{l — z)), yields successively 



E(p) = -ir 



jd'q (l-«/d. 
= -m^ Jd\ {l-tD jdz 



z{l — 2z)m^{p^ — 2pq) 



\w? — z{l — z)q^][rn? — z{l — z){p — qY\[q^ — m 
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zm?' {I — 2z){p^ — 2pq) 



z'^{l — zY — m'^] [{p — qY — rn^] [q^ — rn^] 

In evaluating the integration over the momentum q one introduces two more Feyn- 
man parameters x and y such as to obtain 

-1 .1 r ^2 



T.{p) = -2m^ jdz jdx \y jd'q ' ^1) _ 

{I - 2z){p^ - 2pq) 
{(1 — X — y){q^ — m?) + x[{p — q^ — m?] + y{q^ — m?)Y 
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Translation of the integration variable q ^ q + xp decouples the remaining internal 
momentum from the external momentum p so that one obtains 



rl pi pl—x f- 

E(p) = -2i7r2 dz dx dy d\ (l - ij) 
JO Jo Jo J 

(1 - 2z) (p2 - 2pq - 2xp^) 

z{l — z^ {g^ — x^p^ — (1 — y)m'^ — ym? + xp'^Y 

rl [-1 j-l-x r-oo 2 2 

= {-2m'^f dz dx dy drr"^ ^ 



^0 ^0 ^0 ^(1 
{l-2z){l-2x) 

+ x'^p'^ + (1 ~ y)m'^ + ym? — xp'^}^ 



TT^ dz dx 
Jo Jo Jo 



1 n-x ^2p2 

z[l — zy 
(1 - 2z){l - 2x)x{x - l)p^ 



[x{x — + (1 ~ y)m? + ym?][{l — y)m'^ + ym^ 

rl pi rl—x -j^ 

— TT^ dz dx dy 



JO 



(1 — y)m'^ + z{l — z)ym?' 

2^4 



z{l — 2z){l — 2x)x{x — l)m p 



x{x — '\)z{l — z)p'^ + (1 — y)rn? + z{l — z)y'm? 
The integration over the parameter y yields 
"1 [l-2z){l-2x)p^ 



= TT*^ I dz dx 
Jo Jo 



{z - - z + z^) 



In 



{—X — (1 — x)z{l — z)) (x(l — x)z{l — z)p'^ — w?) 
—xz{l — x){l — z)p'^ + w? — (1 — — 2; + z'^)m/ 



Finally, one notices that the following integral vanishes, by the antisymmetry of 
the integrand under x < — > (1 — x), 

dz fdx ii^Mi^^in ( <^-^W-^)p'-m' 



JO 



{z - - z + z^) V 
Making use of this fact one obtains the result 

6(27r)'^, , 4 A f\ {l-2z){l-2x)p 



= ^i^S(p) [dz [dx 

9 Jo Jo 



[z - - z + z^) 



f [.X + (1 - x)z{l - z)] 

\— x(l — x^z{\ — z)p^ + (1 ~ x)z{l — z)m? + xm? ^ 

This is the result shown in Eq. (23b). 
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